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We study Markovian stochastic motion on a graph with finite number of nodes and adiabatically 
periodically driven transition rates. We show that, under general conditions, the quantized currents 
that appear at low temperatures are a manifestation of topological invariants in the counting statis- 
tics of currents. This observation provides an approach for classification of topological properties 
of the counting statistics, as well as for extensions of the phenomenon of the robust quantization 
of currents at low temperatures to the properties of the counting statistics which persist to finite 
temperatures. 
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I. INTRODUCTION 

This is the second manuscript of the series devoted to 
quantization of currents in periodically driven stochastic 
systems. In the previous article [1], the main object of our 
study was the average current generated per period tjj 
of the driving protocol. The average current, however, is 
only one of the characteristics of the observable stochastic 
currents in mesoscopic systems. The current, Q, gen- 
erated by a periodic driving protocol, can be treated 
as a random variable. In realistic experiments, meso- 
scopic stochastic fluctuations of currents (fluxes) can be 
observed via accumulation of sufficient statistics so that 
one can study not only the mean and the variance but 
also the values of a large number of higher cumulants, 
associated with the current distribution. 

The theory of the full counting statistics of currents 
in periodically driven stochastic systems has attracted 
lot of attention [2H22]- This interest was motivated by 
experiments with molecular motors |23H25j and meso- 
scopic electric circuits HHS], as well as applications in 
soft matter , nonequilibrium thermodynam- 

ics [m I^M5T] and hfe sciences [12 [35]. Therefore, the 
central object of study in the present article is not the av- 
erage current, but rather the complete counting statistics 
of currents, which is represented by the probability dis- 
tribution function V{Q) of currents Q. Specifically, we 
will concentrate on the Large Deviation (LD) function, 
S{Q), that describes the long-time asymptotic behavior 



of the probability distribution function, given by [331 [31] 

r{Q)^e-^^^Q\ (1) 

where N is the number of repeated full cycles of the driv- 
ing protocol. This function is independent of initial con- 
ditions. The universality property, given by Eq. ([T]), is 
known to become exact in the "thermodynamic limit" , 
iV ^ oo. 

In this manuscript we will explore the property of the 
counting statistics that, to our knowledge, has been never 
considered before in the framework of classical stochastic 
Markov chain models. We will show that certain topo- 
logical invariants, in our case the Chern classes, can be 
associated with counting statistics of currents in periodi- 
cally driven systems. We are only aware of similar obser- 
vations in quantum mechanical systems at zero temper- 
ature [351 136j , which originate from essentially different 
initial physical assumptions. 

The first goal of the present manuscript is to explore 
simple models, in which the full counting statistics can 
be obtained explicitly and the corresponding topological 
properties can be clearly illustrated. Our second goal is 
to demonstrate that integer quantization of currents is 
a special manifestation of topological effects in the full 
counting statistics. This relation is important since in 
the first manuscript of the series [T] , we classified integer 
quantized currents. By relating quantized average cur- 
rents to topological invariants of the full counting statis- 
tics we build an approach for classification of such invari- 
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ants. 

The manuscript is organized as follows. In section [TTl 
we introduce the basic notation, including the definition 
of the generating function and its evolution equation. In 
section [Tnj we consider two simple models, one with in- 
teger, and one with fractional average current quantiza- 
tion. We will study in some detail the low-temperature 
limit for these models. In particular, we will obtain the 
generating functions associated with the current distri- 
butions and demonstrate the topological nature of cer- 
tain ingredients in the expressions by relating them to 
Chern classes. In section |IV| we extend the picture of 
the integer average current quantization to the level of 
the LD function S{Q). This will be done for an arbi- 
trary network, whose geometry is described by a gen- 
eral type of a finite connected non-oriented graph. We 
will show that under the condition of the average cur- 
rent integer quantization, the expression for the generat- 
ing function (for an arbitrary graph geometry) involves 
a set n{s) = {na{s)\a G Xi) of integer-valued topolog- 
ical invariants, identified as the Chern numbers. The 
main result of this section is represented by the claim 
that the LD function in the low temperature adiabatic 
limit S{Q) = S^{Q — n{s)) is determined by an even 
function S^{Q), with an integer-valued shift, determined 
by a topological invariant n{s), identified as a set of 
Chern numbers. This result not only provides a Chern- 
class interpretation of the average current quantization, 
but also establishes a stronger property of the current 
PDF, i.e., that the average current is the only odd cu- 
mulant that survives in the low-temperature adiabatic 
limit. In section [Vj we demonstrate how the robust pa- 
rameter subspace Mxi defined in the previous article, 
arises naturally when the LD function is studied in terms 
of the corresponding generating function, with the fo- 
cus on its global properties. We will then develop topo- 
logical interpretation of the (fractionally) quantized cur- 
rents. The conclusion section summarizes our findings 
and discusses possible future research directions. Two 
appendixes serve the goal of providing some very basic 
information on Chern classes, which is an attempt to 
make the manuscript self-contained. 



II. GENERATING FUNCTIONS AND TWISTED 
MASTER OPERATORS 

In this manuscript we show that integer quantization of 
pumped currents in the low-temperature adiabatic limit 
is a consequence of a more general property of the long- 
time asymptotics of the current probability distribution 
ViQ). 

We start with a case of general (not necessarily pe- 
riodic) pumping. The easiest way to study the current 
probability distribution is to introduce the corresponding 



generating function 

Z{\;t) =Y.'P{q;t) [] A^" 

q aeXi 
= ^7'(q;t)e*2:„ex,x..»^ (2) 

<? 

where x — (Xa|o G ^i) is a multi- variant argument 
of the generating function in additive form, represented 
by a vector spanned on the graph links. We denote by 
A — (Aq,|q; G Xi) its multiplicative counterpart; they are 
naturally related by A = e^^. The space in which argu- 
ment A resides is a multidimensional torus (a cartesian 
product over the graph links a of the circles |Aq| = 1) is 
denoted by T. When a link a with da — {i, j} is not a self 
link, i.e., it connects two different nodes i and j, it is con- 
venient to represent Xa by Xia or Xja, with Xia = -Xja- 
For the multiplicative counterpart we use the notation 
Xia or Xja, with XiaXja = 1. The generating function 
Z is a periodic function with respect to parameters Xja 
with period 2tt. Stated equivalently, Z is well-defined 
(i.e., single- valued) function of A. These properties re- 
flect the discrete nature of the transition events through 
the links. 

Similar to the case of the average current, the gener- 
ating function has a useful operator representation 

Z{\;t)= ^ (tcxp( f dt'Hxit'))) p,(0),(3) 
where the twisted master operator Hx(x) has the form 

m 

+ '^'^\xa + Xa'-2)e^(^^-^^)p,. (4) 

The word "twisted" here reflects the way the operator 
Hx is constructed, namely it is obtained from the mas- 
ter operator by multiplying (twisting) the rates kj^ in 
its off-diagonal components, responsible for the transi- 
tions through the links a, with the factors Xja- The 
diagonal elements of the twisted master operator H\ are 
the same as in iJ, i.e., twisting affects the off-diagonal 
components only. Twisting allows the counting of jumps 
through links. Note that the relations XiaXja = 1 for 
the links a with da = i.e. connecting the nodes i 

and J, results in counting the fluxes, i.e. the differences 
of the numbers of jumps along the links in two different 
directions. Also note that relaxing the above conditions 
would result in a more detailed generating functions that 
allows counting the number of jumps along any link in 
each direction separately. However, in this manuscript 
we focus on the fluxes, and hence the above relations are 
always in place. 
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Direct inspection of the matrix elements {Hx)ij of the 
twisted master operator shows that, provided |Aq| = 1 for 
all links a € Xi, or equivalently X E T, the operator Hx 
is Hcrmitian with respect to the weighted scalar product 
in the population vector space, defined by 

and we reiterate that, according to the notation, intro- 
duced in the first manuscript of the series Kj = e^^^ 
and ga = e~^^°', so that the rates are parametrized by 
kja = ffaKj. The Hermitian nature of the twisted master 
operators Hx is important for application of the adia- 
batic limit, since it implies that all eigenvalues of Hx are 
real so that that the notion of a ground state as an eigen- 
vector with the lowest absolute eigenvalue is in place. It 
is also important to note that the untwisted master op- 
erator H = Hi is also real, and hence its eigenmodes, 
i.e. the Boltzmann distribution and relaxation modes 
are represented by vectors with real components. The 
eigenmodes of a twisted operator Hx in the generic case 
A ^ 1 are represented by vectors with generally complex 
components. It is the latter circumstance that allows 
the topological invariants, namely the first Chern class 
to participate in the game. 

In the case of periodic driving, we choose the observa- 
tion time, to = Ntd, to be an integer multiple N of the 
protocol time period td, and note that for large N the 
expression for the generating function [Eq. ([3|] is domi- 
nated by the eigenmode of the twisted evolution operator 

Ux{s;td) = Texp (^td dTHx{s{T))^ , (6) 

with the largest absolute value of the eigenvalue, the lat- 
ter denoted by zx, so that we have for the generating 
function in the large N limit with the exponential accu- 
racy 

Z(A; A/-) - zj^ = e^"\ (7) 

where we have introduced the cumulant generating func- 
tion ujx defined by zx = e^^ . The LD function, S{Q), 
is related to wx via the inverse Legendre transform, 
^iQ) = maxA {~xQ + ^x)- Since A = e'-*', we may write 
uj-^ instead of ujx with the same meaning. 

The generating function zx possesses a gauge invari- 
ance that reflects conservation of currents at long times. 
A gauge transformation, naturally labeled by a vector 
if = {^Pj\j G ^o) in the additive representation and 

= £ -^^0)7 with hj — e*-^J in its multiplicative 

counterpart, transforms x to x' and A to A', so that for 
a link a with da = we have x'ia = Xia - + Vj, 

and A^^ = h~^Xiahj, respectively. The gauge invariance 

z^> = Zy^, ZX' = Zx (8) 

follows from an obvious fact H-^'^ij = e~'^^' H-^^ijC^^^ , 
or equivalently Hx'.ij = Hx.ijhj which does not 



change the eigenvalues of the evolution operator, given 
by Eq. ^. This can be verified directly, and refiects 
the law of current conservation (or, in other words, the 
continuity equation) 

dQ = 0, (9) 

the latter representing a discrete counterpart of the well 
known equation divQ — 0. 



III. TOY MODELS 

The goal of this section is to illustrate the general con- 
cepts using some simple examples. To that end we con- 
sider in detail two simple models. The first one that 
contains two nodes connected by two links has been con- 
sidered in the first manuscript of the series to illustrate 
the main concepts relevant for the average current quan- 
tization. The second model extends the first one to the 
case of two nodes connected by an arbitrary number N 
of links. In this section we calculate explicitly the gen- 
erating function Z(X; N) and its long-time LD counter- 
part Zx [related by Eq. ([T])] for these simple models and 
demonstrate in this simple case how the average current 
is related to the topological invariants of the full counting 
statistics. The general case is considered in section |IV| 
Since in the physics literature it is common to use the ad- 
ditive representation for the argument of the generating 
function, throughout this section the additive argument 
X will be used, whereas in section |1V| we switch to the 
multiplicative argument, as more appropriate and conve- 
nient for dealing with the case of a general network. We 
reiterate that the multiplicative representation is related 
to its additive counterpart by A = e^'^. 



A. Markov-chain on a 2-nodes-2-links graph 

Our simplest model with a pair of nodes connected by 
a pair of links is model in Fig. [ij a particle randomly 
jumps between two sites along two different paths. The 
vector p — (pi,p2)i of the node populations pi and p2 
evolves according to the master equation, 

with the master operator 

^ f -i^iigi + 92) K2{gi+g2) \ /j^N 
V '^1(51 +52) -^2(91+ 92) J' 

where the kinetic rates are expressed in terms of the node 
energies Ej and barrier heights Wa by kja = Kjga, Kj — 
9a = e~^^°', with /? = l/{kBT) being an inverse 
temperature. 
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FIG. 1. The "2-nodes-2-links" model of stochastic transitions 
between two states along two possible paths. Different paths 
are characterized by the values of the corresponding poten- 
tial barriers Wi and W2- The two states are characterized by 
the values of their well depths, Ei and i?2. Periodic modu- 
lation of these parameters leads to appearance of a nonzero, 
on average, circulating current in the preferred (clockwise or 
counterclockwise) direction. 



The main subject of our studies presented in the first 
manuscript of the series, was the average current per pe- 
riod, given by 

Q = N-'^ / dtJ{t), N ^oo. (12) 
Jo 

In the present manuscript we rather treat both Q and 
J{t) as vectors with fluctuating components. Note that, 
in the toy model, current conservation requires that Qi = 
Q2 = Q, so it is enough to study the current component, 
Q, only through the first link. 



B. Generating functions and twisted master 
operators 

We will study the currents that circulate in a coun- 
terclockwise direction, see Fig [T] under the action of a 
periodic driving of parameters. The currents are stochas- 
tic and, generally, can be characterized by probabilities, 
Pj.n{t), that at time t the system is at node j, having 
performed n transitions during the time segment [0, t] 
through the link 1 (with the barrier Wi ) , counting coun- 
terclock/clockwise transitions with the ± sign, respec- 
tively. 

The Master Equation for Pj.n{t) is given by 



Pl,n = -(fell + fcl2)Pl,« + ^2lP2,(n-l) + fc22P2,n, 
P2,n = -(^21 + fc22)P2,n + + fcl2/'l,n- 



(13) 



Rec allin g that kja = i^-jSa, followed by multiplying 
Eq. ( 13 ) by e'-^" and further performing the summation 



over n we find 



where 



Zj=Zj{x,t)= 



Pi 



(14) 



i = l,2, (15) 



and where 



Ki(5ie'^+52) -1^2(91+92) 



(16) 



where is the twisted master operator, in accordance 
with its general definition, given by Eq. Q. 

After many periods of the driving protocols the infor- 
mation on the initial state of the probability vector is 
dissipating, and the information on the currents per pe- 
riod of the driving protocol is contained in the largest 
eigenvalue of the evolution operator 



[A, = T exp 



dtH^{t) 



(17) 



with the twisted master operator, given by Eq. (16), so 



that after N periods of driving, and for large N , the 
generating function of currents is given by 



z{x-N) ^ z^{x:N) + Z2{x;N) ^ z 

which represents a particular case of Eq. ^. 



(18) 



C. Symmetries of the twisted master operator and 
counting statistics 

The closed nature of the considered stochastic process 
implies that the particles are not entering and leaving 
the network, there is actually one particle, jumping be- 
tween the graph nodes along the links, imposing cur- 
rent conservation in the long-time limit. The Arrhenius 
parametrization of the transition rates kja that reflects 
"instantaneous" detailed balance implies that the only 
non-equilibrium feature of the process that leads to the 
current generation is the time-dependence of the rates. 
The two aforementioned features of the stochastic pro- 
cesses, studied in this series, result in the corresponding 
symmetries of the twisted master operators H-^, more 
specifically, gauge invariance and Hermitian nature, re- 
spectively. These two symmetries of H-^, which have im- 
portant implications for the properties of the generating 
function, were already outlined in section [IT] for the gen- 
eral case of an arbitrary network. In this subsection we 
discuss the gauge invariance and Hermitian of the twisted 
master operators in more detail and more explicitly for 
the case of our simple toy model. 



1. Current conservation and gauge invariance 

The node populations cannot increase with time indefi- 
nitely, hence the currents, generated over a large number 
of cycles of the driving protocol, must satisfy the con- 
dition that the difference of such currents entering and 
leaving each of our two nodes is zero. If the currents 
through both links of our two-node model were indepen- 
dent, we should have to introduce a two-component vec- 
tor X = (xi,X2) of the counting parameters and a more 
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general twisted master operator 



-K-iigi + 92) 



^2(316" 



Kl (516*^1 + .926 -^2(31+52) 



(19) 

compared to the one given by Eq. (16). However, due 



to current conservation, the evolution with the opera- 



tor given by Eq. ( 19 ) contains the same information in 



the long-time limit as the simpler evolution operator of 
[Eq. ( 16 1. This can be justified by the fact that the eigen- 



value z-)^ of the evolution operator Hj^ is gauge invariant, 
i.e., it stays unchanged under the gauge transformations. 
A gauge transformation Xi — Xi ~ 'fii + ^2 and simulta- 
neously X2 — X2 — ^2 + ^1, parametrized by the variables 
Lpi and <y92, leaves z-^ invariant, i.e. 



which follows from the fact that 



with 



e^'^i 
e*' 



(20) 



(21) 



(22) 



The transformation, given by Eq. (21), does not change 



the eigenvalues of the evolution operator [Eq. (17)], since 
if \u) is an eigenvector of the evolution operator with the 



twisted master operator Hy [Eq. (21)], then 



\U) IS 



191] 



an eigenvector of the evolution operator H-^/ [Eq. ( 
with the same eigenvalue. This gauge symmetry al 
the parameters ipi and ip2 to be chosen in a way that 
fixes the gauge with x'2 = 0- 



2. Detailed balance and Hermttian property of twisted 
master operators 

When detailed balance is imposed, i.e., the rates 
satisfy the Arrhenius parametrization and are time- 
independent, no average currents are generated in a 
steady state. In fact, all odd cumulants of the current 
distribution must be also zero due to the time-reversal 
symmetry, which is guaranteed by the detailed balance 
conditions. 

In the regime, in which all odd cumulants are zero, the 
generating function should be purely real, meaning that 
the eigenvalues of the twisted master operators with ki- 
netic rates constrained by detailed balance must be also 
real. In fact, one can check directly that both eigenval- 
ues of the operator, given by Eq. (16), are real- valued 
functions of %. 



At this point, we recall that the eigenvalues of Hermi- 
tian operators, which play a crucial role in quantum me- 
chanics, are all real. Apparently, the matrix in Eq. ( 16 ) 



does not coincide with its complex conjugate. Neverthe- 
less, a simple transformation [34l 



I3E/2 



(23) 



makes it Hermitian, and if |u) is the eigenvector of the 
Hermitian matrix H!^ [Eq. (23)], then e"'^^/2|u) is the 



eigenvector of the twisted master operator [Eq. (16)], 
and the eigenvalues of both operators are identicaL In 
other words, one can say that the twisted master operator 
is Hermitian with respect to a scalar product, weighted 
by set of additional energy-dependent factors [Eq. ([5|]. 



D. Topological invariants of full counting statistics 

Depending on the values of the parameters, the twisted 
master operator in Eq. ( 16 ) may have degenerate eigen- 



values. It is a straightforward exercise to show that eigen- 
value degeneracy happens when three conditions are si- 
multaneously satisfied |15j . namely, 



El = E2, Wi = W2, X 



(24) 



Applying the language developed in the first 
manuscript of the series we say that parameters are bad if 
they satisfy conditions Ei = E2 together with Wi = W2 ■ 
Otherwise we will say that parameters are good. We recall 
that the space of good parameters is denoted M = M\D, 
where the discriminant space D is the bad parameters. 
For our simple toy model the space — Ai \ of ro- 
bust parameters, also introduced in the first manuscript 
of the series coincides with the space of good parame- 
ters. Since the good parameter space M and its discrimi- 
nant subspace D are 4- and 2-dimensional, respectively, a 
generic driving protocol avoids the bad parameter space, 
and generally one needs certain symmetry in the system 
for a driving protocol to go through the double degener- 
acy point, other than by just an accident. 

In this subsection, we show that the full counting 
statistics of currents in our 2-node-2-link model, driven 
along a closed path in the good parameter space, can be 
classified using an integer number. 

Since the twisted master operator is periodic in x, 
we can view x as a parameter that resides in a cycle 
[0,27r]. Due to the periodic nature of driving, we can 
also think about the dimensionless time, t, as a variable 
that parameterizes another cycle [0, 1]. Consider a torus, 
— X S^, parametrized by {x,t) with the points 
X = and x = 27r, as well as r = and t — 1 identified. 
With each point on the torus, we can associate an eigen- 
vector, |u^(t)), of the operator that corresponds to 
the larger eigenvalue (''')' i-^-' 



^xMI%M)=exMI%W)- 



(25) 



does not look Hermitian in the sense that its transpose 



The eigenstate choice is not unique because the multi- 
plication of the eigenvector by a complex nonzero factor 
Cxi''') provides a valid eigenvector. 
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Consider a space, referred to as a fiber bundle, that 
is made of the base space represented by the torus, 
T2 = 51 X 5\ with each point (x, r) e of the base 
equipped with a fiber, represented by complex nonzero 
factors Cx(''') ^^at describe the freedom in the choice of 
an eigenvector. During adiabatic evolution of param- 
eters, the twisted population vector that satisfies the 
twisted mater equation [i.e., the one obtained by replac- 
ing the master operator H{t) with its twisted counter- 
part H^{t)] follows the path \u^{t)) of the instantaneous 
eigenstates at the instant values of (x, t). Since we chose 
the contour to be always in the space of good parame- 
ters, this evolution leads to a unique continuous evolution 
of factors Cx ('''); i-^-j fo'^ ^ fixed x, this evolution corre- 
sponds to some curve in the fiber bundle. The integer 
nature of the number of jumps requires the argument 
X of the generating function to be a periodic variable, 
in particular Cx+ztt = Cx- known that fiber bun- 

dles with continuous fiber fields Cx (''")) taking values in 
nonzero complex numbers, are characterized by the first 
Chern class, and in the case of the base being a com- 
pact oriented 2-dimensional manifold, e.g., a torus 
the first Chern class is represented by an integer num- 
ber, hereafter referred to as the Chern number. 

In our case the Chern number has a simple and trans- 
parent interpretation. After a large number of driving 
cycles the information on the initial state of our stochas- 
tic system is lost and the vector of generating functions 
satisfies the relation, \u^{t + 1)) = z^\u^{t)), where we 
denote = Cx(l) 7^ ^- ^'^^ periodicity in x, we 

can classify this function according to how many times 
the curve Ztj. winds around the zero point in the complex 
plane while parameter x grows from zero to 27r. The 
Chern number in our model corresponds merely to this 
winding number. 

Having established that the counting statistics of cur- 
rents in periodically driven systems can be classified by 
topological invariants, the next questions are: 

(i) How to calculate such invariants without solving a 
time-dependent twisted master equation explicitly? 

(ii) What are the measurable manifestation of topolog- 
ical invariants in stochastic processes? 

We will address these two questions, for the toy model, 
in the following two subsections. 



E. Derivation of the winding number 

In this subsection we establish a connection between 
the geometric phase of the generating function for our 
toy model and the winding number, and compute the 
latter using its topological nature by applying the low- 
temperature limit. 



1. Geometric phase and winding number 

In [TOj [13] it was shown that in the adiabatic limit the 
cumulant generating function (CGF) uj^ can be written 
as a sum of two contributions, 

^x=<+^^- ^ = 44 (26) 

referred to as the geometric and quasisteady terms. The 
quasisteady contribution has a natural form 

w«= / dte^{t)=TD dre^ir). (27) 

The geometric contribution depends only on the ge- 
ometry of the contour in the parameter space M. (which 
justifies the notation), and is independent of the rate of 
motion along this contour, as long as the adiabatic ap- 
proximation can be employed. It has the form 

^x = - / Mx,x) ■ dx, A{x,x) = {u^{x)\d^Ux{x)), 

J S 

(28) 

or equivalently 

^X^- f dTAo{x,T), Ao(x,t) = {u^{T)\drU^{T)) . 

Jo 

(29) 

Introducing also 

A,{x,T) = {u^iT)\d^u^{T)), (30) 

we can interpret A — {Aq,Ai) as a vector potential 
(or equivalently an abelian gauge field) defined on our 
torus. The reason why A is referred to as a gauge field 
originates from the fact that the ground state \u{x,t)) 
is defined up to a scalar factor Cx (''')' ^^'^ ^ rescaling 

= Cx(^)l«x(^)) = e^^^'"Vx(^)> leads to a gauge 
transformation of the vector potential 

A'j = Aj - djcl), \u')^e^\u) (31) 

where do = and di = d-^. 

In terms of the cumulant generating function (CGF), 
the topological invariant has a new physical interpreta- 
tion: while the generating function z-^ is a periodic func- 
tion of Xi the CGF, w^, may not be necessarily peri- 
odic, but it rather has the property, ij-'x+27r = + 2'KiQ, 
where Q is the integer that coincides with the winding 
number. The CGF is a measurable function. By measur- 
ing sufficiently many cumulants of the distribution, one 
can determine Q if the corresponding representation of 
the CGF in terms of a series of cumulants converges for 
arbitrary x € [0, 27r]. 

The geometric contribution does not appear in steady 
states, since it is nonzero only when the parameters are 
varying along the contour s. The steady state CGF 
would be given by s^td, which means that is just 
the time average steady state contribution. The eigen- 
values of the twisted master equation are real because of 
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the Hermitian property of the twisted master operator 
at detailed balance conditions on the kinetic rates. This 
means that the strict quasisteady state contribution to 
the CGF is purely real. In particular, it contributes only 
to the even cumulants of the currents and is not related 
to the winding number. Thus the information about the 
topological invariant is contained in the geometric con- 
tribution, w^, specifically, the winding number, Q, given 
by 



Q 



1 

27ri 



1 

2TTi 



dx 



8 z3 



(32) 



Note that since is a single-valued function of x, one 



can replace oj^ and z| in Eq. ( 32 1 with uj^ and Zy, , re- 
spectively. 

There is a delicate and important detail that needs 
to be discussed. The ground state |u^(t)) is globally 
defined on the torus only up to a scalar position- 
dependent factor Cxi"^) ^^^^ may or may not be fixed 
globally in a smooth (or even continuous) way. Actually, 
the nonzero value of the winding number Q forms a topo- 
logical obstruction to fixing Cx(''') globally. Therefore, as 
explained in some detail in appendixjX] the vector poten- 
tial A = (Aq, Ai) is defined globally only up to a gauge 
transformation [Eq. (31 1]. More precisely this means that 



the vector potential can be defined in some regions that 
cover our torus, so that in the intersection of any two re- 
gions the vector potentials differ by a gauge transforma- 



tion. Therefore, the expression for a;| [Eq. (29)] should 



be understood in the following way. Fixing x we can 
define \uy^{j) together with the scalar factors on a cycle 



< T < 1 and further apply Eq. ( 29 1 . Being represented 



by an integral of the vector potential over a closed con- 
tour, is gauge invariant and therefore well-defined. 

It is still possible to come up with an integral represen- 
tation of the winding number Q with a gauge-invariant 
integrand. To that end we introduce the curvature 
( "magnetic field" ) F of the vector potential 

F(x,r) = 5oAi - ai^o = (dru\dxu) - (33) 

sometimes referred to as the Berry curvature |13j . so that, 
as justified with some detail in appendix [A| Eq. (32) can 
be recast as 



1 

2ti 







2ir r-\ 

dx / dTF{x,T), 
Jo 



(34) 



with the curvature, due to its gauge-invariant nature, 
being a smooth function globally defin ed on our torus 
T^. Also note that the r.h.s. of Eq. (34) actually re- 



produces the original definition of the first Chern class 
for 2-dimensional manifolds, so that Eq. (34) provides a 



Chern number view of the topological invariant Q of the 
full counting statistics. 



A calculation of our topological invariant via Eq. ( 34 ) 



can still be technically a hard problem. It requires know- 
ing the eigenstates of the twisted master operators, that 



might not typically be found analytically. Fortunately, 
for stochastic models there is an alternative procedure. 

A key observation that makes calculation of the wind- 
ing number in our toy model easy is that the degeneracy 
conditions [Eq. (24|] do not depend on the temperature, 
so if we consider a continuous class of periodically driven 
systems that differ only by temperature, the winding 
numbers that characterize each of these models should 
be the same because continuous changes of temperature 
can lead only to continuous changes of Cxi''')^ ^^^^ cannot 
change topological invariants of the fiber bundle. Thus 
we can, in fact, calculate the winding number in the limit 
of zero temperature, /? — ^ cx3. 



2. Winding number in the low-temperature limit 

Taking the low-temperature limit leads to considerable 
simplifications since the ground state of the master oper- 
ator can be found explicitly at each point (x, r) € of 
the torus up to exponentially suppressed corrections. For 
our two-state model, whenever the barriers, Wi and W2 
are different, we can disregard the factors which contain 
e-PWi or e-^^\ if Wi > W2 or W2 > Wi, respectively 
Still, there can be short time intervals when the driving 
protocol passes through the point Wi = W2 ■ Precisely at 
the latter intervals, the restriction of dealing only with 
those protocols that lie entirely in the good parameter 
space leads to the condition Ei ^ E2, and hence we can 
disregard terms containing e^^'^ or e'^^^, if, respectively. 
El < E2 or E2 < El. 

Consider again the driving protocol, E2 = W2 — 0, 
Ei{t) — cos(27rr), Wi{t) — sin(27rr), that has been al- 
ready explored in the first manuscript of the series to 
study average currents. One can distinguish four time 
segments with distinct behavior of the system: 

(a) for T e [— e, e], where e < tt/2, the path goes through 
the point Wi = W2 but we always have Ei > E2; 

(b) for r e (e, TT — e), we have Wi > W2; 

(c) for T £ [tt — £, TT -|- e], we pass through the point 
Wi — W2 but we always have Ei < E2; 

(d) for r e [tt + £, 27r — e], we have W2 > Wi. 

Consider stochastic dynamics on each of those seg- 
ments in the low-temperature /? — > 00 limit. At segment 
(a), we can disregard the terms with the factor e^^^, so 
that 



X 



-Ki{gi+g2) 
Ki(5ie*x -1-52) 



has a right eigenvector 



(35) 



(36) 



that has a simple interpretation: the fact that it does not 
depend on time means that it represents a process during 
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W. 



El 



FIG. 2. Two steps (a) Wi < Wi and (b) VFi > VK2. At the 
low temperature limit, the effective model is a diffusion on a 
tree graph that consist only of one link and two sites. 



which the system remains in the second node without 
producing any current. 

By analogy with (a), the evolution at segment (c) is 
also trivial, and it is represented by the eigenstate 



(37) 



indicating that the system stays at the first node for the 
entire time segment. 

At segment (b), we can disregard the terms with the 
factor e~^^^ . The approximate twisted master operator, 
for W2 < Wi, turns out to be independent also of the 
counting parameter 



X 



„/3 (£2-1^2) 
_„I3(E2-W2) 



(38) 



It is equivalent to the non-twisted master operator of the 
model in Fig. [sjb). The largest eigenvalue of such an 
operator is zero. The right and left eigenstates of H^^ 
that correspond to the zero eigenvalue are given by 



--PE2 



,(6)1 = 



(1,1), 



(39) 



with Z = e^^^^ + e~^^^ being the the normalizing fac- 
tor. Generally, the eigenstates of the operator given by 
Eq. (38) can be chosen with an arbitrary phase e"^^-^'^^. 
However, in order to apply Eq. (29 1, we should choose 



the eigenstates of the twisted master operator so that 
they continuously change along the contour s. The choice 
of the real prefactor in Eq. (39) guarantees that at the 
points where the segments (cijand (c) overlap with (b), 
the eigenstates coincide up to exponentially suppressed 
terms which may be disregarded. 

Finally, at segment (d) the terms that contain 
can be disregarded, so that 



-I3W2 



H(d) 

X 



_^P(Ei-Wi) 
^PiEi-Wi)+ix 



^l3(E2-Wi)-ix 
„g/3(B2-Wi) 



(40) 



The matrix in Eq. ( 40 ) corresponds to the twisted mas- 



ter operator for the Markov chain model illustrated in 
Fig. |2ja) that has two sites and one link only. Such a 
matrix should have the lowest eigenvalue independent of 
X, because otherwise the tree graph would be able to 
support a steady state current with growing cumulants. 



which is impossible. Since x = corresponds to no twist- 
ing, the aforementioned largest eigenvalue should be ex- 
actly zero. The above intuitive arguments can be put 
on firm ground by noticing that the twisted operator in 
Eq. ( 40 ) can be transformed to an untwisted counterpart 

The 



r(d) 



by a gauge transformation in a form of Eq. (|2l|) 
eigenstates of 
are 



corresponding to the zero eigenvalue 



o-PEi 
,-0E2+ix 



i4>( 



(41) 

The phase = '/'(x, ''") should be chosen in a way such 
that \u^(t)) is continuous at r = 27r — e and t — tt + e 
that corresponds to the intersection points of segment (d) 
with the segments (a) and (c), respectively. This requires 
(f){x, 2tt — e) = and (f){x, n + e) — — x, respectively. 

In such a gauge, the eigenvectors smoothly change from 
one interval into another. Calculating the relevant com- 
ponent Aq of the vector potential, we find that in our 
gauge Ao{t) = at the segments (a), (b), and (c). At 
segment (d), we find 



Aq — idr4>- 

Performing integration along segment (d), we find 



(42) 



drAo(r) 



= -^(0(x,27r-e)-(/.(x,7^ + £)) =zx- (43) 

We conclude that in the low temperature limit the ge- 
ometric phase contributes to the first cumulant of the 
current statistics only, provided that the exponentially 
suppressed corrections are neglected. Another observa- 
tion is that a;| is not a periodic function of x, and the 
winding number is equal to one: 



= 1. 



(44) 



The above arguments can be extended to a general 
periodic protocol in a straightforward way. For our toy 
model, we can always split any periodic driving protocol 
into alternating time segments such that (a) Wi « W2 
while El < E2, (b) Wi>W2, (c) Wi « W2 while E2 < 
El, and (d) Wi < W2- The continuous choice of the 
eigenstates in each of those segments is the same as in 
the above example. Then whenever we find the path in 
the parameter space that goes from (c) to (a) via (d), the 
geometric contribution to CGF acquires a term ix and 
whenever the path goes from (a) to (c) via (d) the CGF 
acquires a term —ix- So the CGF per period of a driving 
is given by 

^x = iQx + / dte^{t) = iQx + td dTe^{T){'i5) 
Jo Jo 

where Q is the integer number, which is equal to the 
number of transitions through the interval (d) taken with 
a proper sign as discussed above. 
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We already mentioned that e{x, t) is purely real and 
contributes to the even cumulants only. The coefficient 
of x" hi the Taylor expansion of the CGF corresponds 
to n-th cumulant of the corresponding probability distri- 
bution. Our results show that in the /3 — >■ oo limit the 
geometric contribution is linear in x- Therefore, in addi- 
tion to identifying the first cumulant, i.e., average current 



with the Chern number Q, Eq. (45) also demonstrates 



vanishing of all higher odd cumulants of the generated 
current in the low-temperature limit. 

In the following sections, we will show that the quan- 
tized current in a generic network is also a manifestation 
of topological properties of the full counting statistics of 
currents in periodically driven networks. 

The main difference between the toy model and a 
generic network is that the latter typically has more than 
one independent current, i.e., the argument of the gen- 
erating function is typically multi-variant. Also, for a 
generic graph, the condition for the largest eigenvalue de- 
generacy is typically temperature dependent. However, 
at sufficiently low, yet finite temperature, the degener- 
acy can be avoided for all values of the CGF argument 
X, provided the parameters belong to the good parame- 
ter space. Indeed, by restricting the twisted master op- 
erator to the minimal spanning tree, which can be jus- 
tified by the same arguments as for steps (b) and (d) 
in the toy model, we find that the largest eigenvalue is 
zero and non-degenerate. Therefore, if the temperature 



is finite, but small enough, the largest eigenvalue is still 
non-degenerate, provided the parameters are good, and 
the topological invariants, obtained by considering de- 
terministic evolution in the context of average currents, 
characterize the full counting statistics, at least in some 
range of temperatures. 



F. Toy model for fractional quantization 

According to [H] , when permanent degeneracies of the 
node energies and/or barrier heights take place, the cur- 
rent, generated by adiabatic periodic evolution of the 
system parameters, is generally fractional, rather than 
integer- valued in the low-temperature limit. In this sub- 
section, we consider a simple model that only slightly 
generalizes the model, considered so far in this section 
(see Fig. [T]), by replacing each of two links by a group 
of links with the permanently degenerate barrier heights, 
as illustrated in Fig. [s] Let = e~^^^ for a = 1, ... M, 
and Qa = e^^^^ for i = Af 1, . . . TV. We also introduce 
a vector % = {xi, ■ ■ - Xn) of counting parameters that 
monitor the number of counterclockwise minus clockwise 
transitions through each link. The evolution operator for 
our A^-component moment generating function is given 
by 



N 

-Ki E 9i 



M 

1=1 



M 



N 



N 

j=M+l 
N 

K2 E 9i 
1=1 



(46) 



r 



In the adiabatic limit, the calculation of the CGF is 
reduced to the eigenvalues and eigenvectors of the 2x2 



matrix in Eq. ( 46 1 , which can be completed explic- 



itly. Moreover, in the low-temperature limit, one can 
split the driving protocol period time into time segments 
intervals (a)-(d), similar to what has been done for the 
2-state-2-link model, and apply similar approximations, 
based on disregarding the links with the higher barriers, 
when computing the generating function, as illustrated 
in Fig. |4] For the same driving protocol, as considered 
for the 2-state-2-link model earlier in this section, the 
calculation, described above, results in 



'E 



M 



M 



'E^Lm+i^' 
N - M 



(47) 

where e(x,T-c) is a non-universal quasisteady-state con- 
tribution that is purely real with e(x,T£)) = e(— X,T£)), 
hence obviously contributes to the even cumulants only, 
and, in this sense, is indistinguishable from the thermal 



noise at equilibrium. The most interesting feature of the 
counting statistics arises from the pre-exponential factor 
that describes the effect of the geometric phase. When 
the driving protocol encloses the degeneracy point in the 
parameter space Q times, the geometric factor in z-^ be- 
comes 

/v^A/ i.,.\Q /^N i^..\Q 

E,=i ' 



E 



M 



N 

j=M+l 

N-M 



(48) 



We reiterate that the contributions to the odd cumulants 
are provided by the geometric factor [Eq. pi] )] only. The 
components of the average generated currents can be ob- 
tained by taking the derivatives of In (z^(r£i)) with re- 
spect to the corresponding counting parameters ixj, us- 
ing Eq. (48), followed by setting x = to zero: 



{Q^)=Q/M, if t = l,...M, 



(49) 



and 



{Q^) ^Q/iN - M), ff i^M + l,...N. (50) 
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FIG. 3. A graph with 2 nodes and N Unks. Nodes are char- 
acterized by energies Ei and i?2. There are two sets of hnks 
that are permanently degenerate, so that links with numbers 
1, . . . Af have the same barrier size W2, and the links with 
numbers M + 1, . . . A'^ have barrier size Wi. 




(a) 



(b) 



FIG. 4. Effective kinetic models that describe evolution of 
the system in Fig. |3] in the limit of low temperature for (a) 
the case of Wi < W2 and (b) the case of Wi > W2 • 



As expected, at low temperatures the average currents 
through individual links are fractionally quantized. The 
fractional quantization is expressed through the integer 
parameters, Q, N and M. The latter two originate from 
degeneracy numbers for the two groups of links. The 
integer number Q is the winding number, which still has 
a Chern class interpretation. This can be seen from the 
fact that if we set xi = ■ ■ ■ = Xm = and xm+i = . . . = 
Xn = X the twisted master operator in Eq. ( 46 1 is the 



same as in the 2-state-2-link model, in which Q has a 
Chern number interpretation. 

We conclude this section with noting that we explic- 
itly calculated the CGF or related quantity z-^ — e^^ 
in the low-temperature adiabatic limit for two toy mod- 
els, both containing two nodes. The result for the 2- 
link model that supports integer quantization [Eq. (45 1] 
looks universal, with the geometric contribution totally 
described by the topological invariant Q represented by 
the first Chern class. Eq. (45) suggests that the above 
universality should hold for a general network. Indeed, 
we will demonstrate in section IIVI that this is the case. 
Therefore, the Chern class interpretation of the average 
integer-quantized current is a consequence of a funda- 
mental property of the CFG in the case of no perma- 
nent degeneracy. The expression for the geometric factor 
[Eq. ( 48 )] in the CFG of the extended toy model that sup- 
ports fractional quantization is more complicated. There- 
fore, although in section|V]we provide a Chern class inter- 
pretation of fractionally-quantized current, based on the 
averaging concept, formulated in the first manuscript of 
the series, the presented interpretation, as opposed to the 



integer case, is somewhat artificial and does not reflect 
any universal topological properties of the CGF, which 
are yet to be determined. 



IV. GENERAL NETWORKS 

In this section, we study the generating function on 
arbitrary graph with the focus on the low-temperature 
adiabatic limit. 



A. Generating function in the adiabatic limit and 
topological contributions 

We start with considering the adiabatic r/j — > 00 limit 
followed by applying the low-temperature /3 — 00 limit. 
In what follows, we will be using the multiplicative repre- 
sentation A for the multi-variant argument of the gener- 
ating function. We reiterate that the components in the 
multiplicative and additive representations are related by 

Adiabatic calculation of the eigenvalue za(s;td,/3) of 
the time-ordered exponential operator Ux{s;tz)^ (3) de- 
fined by Eq. (|6| becomes standard when the ground 
state of the twisted master operator Hx{s;f3) is non- 
degenerate at all times. We reiterate that, as justified 
in section [nj the notion of a ground state is well-defined 
due to Hermitian property of the twisted master oper- 
ators Hx, which restricts its eigenvalues to be real. To 
reproduce the complete generating function in this stan- 
dard manner, we want this non-degeneracy to take place 
for values of the argument A. Stated differently, we would 
like our driving protocol to stay inside the parameter 
subspace MxW) C Mx (that naturally depends on /3) 
characterized by the property that Hx{x; (5) has a non- 
degenerate ground state for all A, as long as a; G M^xiP)- 
The subspace M.^x{P) '^iH be referred to as the non- 
degenerate ground-state space, associated with /3. The 
parameter sets x g M\[j3) that belong to this space and 
the driving protocols that lie entirely inside this space will 
be referred to as non-degenerate ground-state parame- 
ters and driving protocols, respectively. Although at this 
point it would seem that the subspace M.'^xiP) has been 
introduced based on our ability to perform adiabatic cal- 
culation, we will show later that in the low-temperature 
limit the above subspace actually describes the robust 
driving protocols. 

As justified above, for a non-degenerate ground-state 
driving protocol, by applying a standard adiabatic tech- 
nique, we have in the tjj — > 00 limit 



ZA(s;Tz,,/3)«CA(s;/3)e^°^^(«^«, 



A^(s;/3) 



dTUJx{s{T)\(3), 



(51) 



and we reiterate that ujx{x]P) denotes the lowest eigen- 
value of the twisted master operator Hx{x;(3), and 



11 



Cx{s; (3) is the geometric factor. An attempt to repre- 
sent the geometric contribution in exponential form (sim- 
ilar to the second factor) faces a topological obstruction, 
represented by a set n(s; /?) = {na{s; (3)\a S Xi) of inte- 
gers, where na(s; /?) is the index that measures how many 
times (\{s; /3) winds around the zero point G C, while 
Aq, goes over a unit circle |Aa| = 1 (as explained in some 
detail in appendix |A| , whereas all other components of 
A stay unchanged. However, we can still represent the 
geometric contribution 



CA(s;/3) = Cj(«;/3)e^"(^^« 

as a product of a purely topological factor 

(lis; (3)= n A-(-« 



(52) 



(53) 



aeXi 



and an exponential term, where A^(s;/3) describes a 
purely geometric correction to the topological term. 
Summarizing, we have 



zx{s;td,I3) 



aeXi 



and, due to gauge invariance 



da={i,j} 

n,„(s;/3)=0. 



,(54) 



(55) 



which means that the set n{s; /3) of integers represents a 
conserving integer- valued current. 



The obtained adiabatic expression [Eq. (54)] has the 
following important implication on the LD function. Let 
S{Q) be the LD function that corresponds via the Leg- 
endre transform to 



Aa(5;to,/?) = A^(s;/3)+tz5A^(s;/3). (56) 



Then 



S{Q)^S{Q-n{s-/3)), 



(57) 



i.e., the topological factor results in just a shift of the LD 
function and, therefore, affects the average current only. 

We are now in a position to show that the integers 
na{s; f3) that describe the shift in Eq. (57) can be viewed 
as Chern numbers, and, therefore, calculated using the 
integral representation, given by Eq. (A6). Indeed, by 
the definition of the robust parameter space A^x, 
family Hx(x]l3) of twisted master operators provides 
a linear fiber bundle over M.^^^ x T- We reiterate 
that a linear fiber bundle is just a collection of fibers 



parametrized by the points (a;. A) G 



r, 

with the fiber £^ being a complex 1-dimensional vec- 
tor space represented by the ground states \u) that sat- 
isfy Hx{x; (3)\u) = ujx{x] (3)\u) , with ujx{x]j3) being the 
lowest eigenvalue. 

We further observe that a periodic driving protocol s 
defines a set /(s) = {Ja{s)\a € Xi) of 2-dimensional cy- 
cles in M\ X T, associated with the network links. A 



cycle that maps a torus to A^5s: associates with a 
point (r, x) in the torus (we reiterate that < r < 1 
and < X < Stt with periodic boundary conditions) 
the point (s(t), A"(x)), where the components of A"(x) 
are: A2(x) = e*^, and A"(x) = 1 for 7 ^ a. Stated 
less formally, if we view a torus as a product of 
two circles, then the two circles of the cycle fa{s) are 
represented by the driving protocol and the unit circle 
|Aq| = 1, respectively. Restricting the argument A of the 
generating function to the circle A-y = e*-^''"' we identify 
na{s](3) — n\C,{x))] wit h the latter given by Eq. (A8). 
Finally, due to Eq. (|A9| 



n„(s;/3) =ci(£'^)/„(s). 



(58) 



which can be interpreted as that na{s] 13) is a Chern num- 
ber, obtained by evaluating the first Chern class ci(£^) 
of at a 2-dimensional cycle fa{s). In particular it 
can be c alculated by using the integral representation of 
Eq. (lA6). 



B. Low-temperature limit for generating function 

With our focus on the low-temperature limit, it is 
desirable to replace the ground-state non-degenerate 
temperature-dependent subspaces M'^{l3) with some 
limiting subspace ^Ax■ The latter can be naturally de- 
fined by the following characterizing property: x € 
if for some f3o{x) the operator Hx{x; (3) has a non- 
degenerate ground state for all A and all (3 > j3o. The 
subspace M.'^x '^^ ^ very convenient tool for applying the 
low-temperature limit. In particular, for large /? forms 
a family of bundles over A^^f ^ ^ ■ Therefore, they are 
all equivalent (isomorphic), and have the same Chern 
class, i.e., with a minimal abuse of notation we can write 
ci(£). More formally, given a driving protocol s, for 
/3 > max(/3o(s(r))), the Chern numbers Ci(£^)/a(s) are 
/3-independent. Therefore, there is a well-defined low- 
temperature limit 



na{s) = lim na{s\ (3) 

/3— foo 



ci{C)fa{s), 



(59) 



where the rhs means that ci(£)/q(s) = ci[C^)fa{s), cal- 
culated for large enough, still finite (3. 

Integer quantization in the low-temperature adiabatic 
limit can be formulated as a stronger statement in terms 
of the LD function S{Q), or equivalently in terms of zx- 
It includes three statements, (i) The topological term 
Ca(s; /?) has a well-defined low-temperature f3 00 limit 

Cl(^)= n ^""^'^ n„(s)=ci(£)/4s), (60) 

expressed in terms of Chern numbers. This has been 
demonstrated earlier in this subsection, (ii) The purely 
adiabatic contribution Ax{s;f3) possesses time-reversal 
symmetry A_x{s; (3) ~ Aa(s; /3), which is inherited from 
the time-reversal symmetry uj^x{x\l3) = lux{x]/3). The 
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latter property can be validated in various ways. We 
will restrict ourselves to mentioning that physically this 
is equivalent to the even character of the LD function of 
a non-driven system, i.e., when the parameters {E,W) 
do not change with time, (iii) The purely geometric cor- 
rection to the topological term vanishes A^(s; /3) — > in 
the low-temperature /? — oo limit, or equivalently 



hm Ca(s;/3)-C!(s) 



(61) 



Properties (i)-(iii) result in the following form of the 
LD function in the low-temperature adiabatic limit 



S^{-J-P)^S^ (J;/3) 



where {Q; td,/3) is obtained as a Legendre transform 
of Ai{s-J 



In particular, Eq. (62) identifies the low- 

n(s), thus 



temperature adiabatic average current Q{s 
providing a Chern class representation for Q(s 



We further partition the (proper) time segment [0, 1] 
into alternating 0- and 1-segments la = [t^_]^,To] and 
la = [''a:''"a]j respectively, exactly in the same way as it 
was done for average currents in previous article of the se- 
ries. We reiterate that by construction s(Ta), s(T'a) € Uq, 
with Uq and Ui being the parameter subspaces with no 
degeneracy in the energy and barrier data, respectively 
(by definition Aix — UqUUi). So we may introduce the 
nodes ka and k'^, with — ^a, that correspond to the 
minimal values of Ei[Ta) and Ei{T^), respectively. A sim- 
ple analysis shows that on a 0-segment la, in the /3 — >■ oo, 
limit the ground state of the twisted master operator is 
the same as in the untwisted case, i.e., it is localized in the 
lowest-energy node Uj{T) — CaiT)Sj°' , whereas, the sec- 



(62) ond condition in Eq. ( 64 ) , as applied to our case, requires 



Cq to be time-independent. Therefore, we can introduce 
the sets of non-zero complex factors C,a and (a, so that 

k k' 

U]{Ta) = CaSj" and Uj{t^) = Ca^j"", respectively. Note 
that C'a-i — Ca- Then, in the low- temperature limit the 
geometric contribution is represented by the product 



C. Low-temperature calculation of the geometric 
contribution 



(65) 



In this subsection, we perform a calculation of the 
geometric factor Ca(s;/3) in th^ low-temperature limit, 
i.e., derive the expression given by Eq. (61). In the 



first manuscript of this series, we identified the subspace 

Mx — -Mx — D'^x of robust parameters by explicitly de- 
scribing the discriminant space of non-robust param- 
eters. In this subsection we accomplish a more moderate 
task by showing that M.x C A^x; that a driving 
protocol that avoids simultaneous degeneracy of the en- 
ergy and barrier data (we reiterate that this is how the 
subspace A^x is defined) is a robust protocol in the sense 
of low-temperature behavior, and derive Eq. (61 ) for such 



protocols. The validity of Eq. (61) for the driving pro- 



tocols that go through M.x will be instrumental for the 
explicit identification of the broader space A^x- 

We fix a driving protocol s, which is restricted io M.x 
at all times, and the argument A G T of the generating 
function, and perform a direct calculation of the geomet- 
ric factor C,x{s\ /3). By definition, the latter is represented 
by a factor in the relation 



Ki))=Ca(s)I"(0)), 



(63) 



where the time-dependent distribution \u{t)), repre- 
sented by u{t) = (uj(r)|j £ Xq), satisfies the conditions 



Hr)\drUiT)) 



^A(s(r);/3)|ii(r)), 
0, 



(64) 



for all < T < 1. Note that Eqs. (64) represent a sys- 



tem of equations for standard adiabatic theory, with the 
second condition ensuring proper treatment of the geo- 
metric phase (known as the Berry phase for evolution 
with time-dependent Schrodinger equation). 



and the only thing we need to do is to identify the ratios 

aca. 

In the low-temperature limit, the ground state \u{t)) 
on a 1-segment can be found asymptotically exactly by 
replacing the original graph with the relevant minimal 
spanning tree X^^. Physically, we are neglecting the 
transition over the links with the highest values of the 
barriers Wa, which can be done, as long as the graph 
stays connected. The ground-state problem on a tree 
can be solved exactly, since twisting can be eliminated 
by a gauge transformation, while the ground state for the 
untwisted case is represented by the Boltzmann distribu- 
tion. Since evolution with a standard (untwisted) mas- 
ter operator conserves the total population, the second 
condition in Eq. ( 64 ) boils down to the requirement for 



the gauge-transformed ground state to be represented by 
the normalized Boltzmann distribution with some time- 
independent multiplicative factor, i.e., in the original rep- 
resentation on a 1-segment la we have 



(66) 



where h — {hj\j E Xi) represents the gauge transfor- 
mation that eliminates twisting. The latter requirement 
means that for any link a on the spanning tree with 
da — {i,j}, we have A.^^ = Xj^ = h^^hj. It is con- 
venient to choose hk^ = 1, so that the set of above 
constraints has a unique solution: hj is given by the 
product of Aq. over all links a that belong to the path 
hajiXwa^) that connects the reference node ka to node j 
through the minimal spanning tree. Recalling the defini- 
tion of the integer- valued current Qij{Xw), associated 
with the path lij{X\\r), followed by substituting the ob- 
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tained value of h into Eq. (p6|) we arrive at 



(67) 



which, due to pf{E{Ta);/3) -> S^^ and pf (£;«); /3) 
5 J for /5 — ^ oo , results in 



C/Ca = n ^' 



(68) 



Substituting the obtained expression into Eq. ( 65 1 and 



recaUing the expression for the average low-temperature 
adiabatic current Q(s), generated by a good pumping 
protocol, we arrive at 



lim Ca(s;/3)= n 



(69) 



The obtained relation [Eq. ( |69| ] demonstrates the van- 
ishing lim^_j.oo A£(s; /3) = of the purely geometric cor- 
rection in the low-temperature limit, or equivalently the 
statement of Eq. (61), and also identifies 



Q{s) = n{s) 



(70) 



the low-temperature adiabatic integer-valued average 
current as the Chern numbers, associated with the cor- 
responding 2-dimensional cycles f{s). 



V. CHERN-CLASS REPRESENTATIONS FOR 
QUANTIZED CURRENTS 

This section is focused on some aspects o f the Chern 
class representation of Q{s). In subsection VA we re- 



visit the integer case and make an explicit connection 
Mx = -^x — D\ between the extended subspace 
introduced in the context of the generating function in 
subsection |IV B| as the limiting space of the temperature- 
dependent non-degenerate ground-state spaces M\{(i) 
with the reduced discriminant set of the non-robust 
parameters, introduced in the previous manuscript of the 
series for a more general rational case. In subsection |VB| 
we discuss a Chern class representation for the rational 



A. Chern class representation for the integer case 
revisited 



In subsection IV B we introduced the space A^V ^ 



M.X that extends the subspace M.x C M.x of good pa- 
rameters. The extended space A^x been identified as 
a result of studying the low-temperature adiabatic limit 
of the generating function for the full counting statis- 
tics in the integer-quantization case, i.e., the pdf of the 
generated current, more precisely the LD function S{Q). 
In the first manuscript of the series, while considering 



a more general rational case, we extended the subspace 
M^x = M.X — Dx by withdrawing certain unnecessary 
cells of Dx which resulted in a reduced discriminant set 
D^x C Dx- In this subsection we demonstrate for the 
integer case that A^x — — D^x- 

To show that M\ C Mx — Dx, we use a purely 
topological argument to demonstrate that if a point 
X = (E.W) G Di^, i.e., belongs to a cell D^^, with 
T>^ 7^ 0, then x cannot belong to Mx- We reiterate that, 
according to the winding-index picture of quantized cur- 
rents, introduced in section VIII of the first manuscript 
of the series for a more general case of rational quanti- 
zation, further applied to the particular case of integer 
quantization, the discriminant space Dx is represented 
by a union of the closed cells D^, the latter obtained ap- 
plying the topological closure to the corresponding open 
cells D^, labeled by w = (jj',aa',p), where j ^ / and 
a ^ a' is a pair of distinct nodes and links, respectively, 
whereas p is an ordering ai < . . . < < {a, a'} < 
ckr+i < . . . < Q?|Xi| of the graph links. The cell Di^ con- 
sists of all parameter sets {E, W) with the constraints 
Ej = Eji < Ei for any i ^ and Wa^ < . . . < Wa,. < 



< Wo, 



Stated equivalently 
the cell D^ describes the parameter sets with simulta- 
neous lowest-energy degeneracy Ej = Eji and barrier 
degeneracy Wa = Wa', whereas we have a specific or- 
dering of the barrier heights, given naturally Wa = Wa'- 
We also reiterate that V^^ = Q{s) represents the low- 
temperature adiabatic current generated over a driving 
protocol s that winds around a point x G Di^ in its small 
neighborhood without intersecting with the cell D^^. As 
demonstrated in section VIII of the first manuscript of 
the series, the current Q{s) does not depend on a partic- 
ular choice of a point x € Di^, as well as a driving pro- 
tocol s, as long as the latter belongs to a small enough 
neighborhood of x and winds around D^^ once. The de- 
scribed property refiects the topological nature of the 
low-temperature adiabatic limit of the generated current. 

To show that Mx C Mx — Dx we use the proof 
by contradiction concept. Suppose not, i.e., x G Mx- 
Since Mx is open (i.e., also contains some neighborhood 
of a;), and since the cells D^ have codimension 2, we 
can consider a 2-dimensional disc A centered in x that 
is contained in the aforementioned neighborhood and in- 
tersects D^ transversally. Let s be a driving protocol 
that goes over its boundary once. Then s is restricted 
to Mx, and A is spanned on s, as well as ^ C Mx- 
By the winding index representation the generated over 
the driving protocol s low-temperature adiabatic current 
Q{s) = A ■ Vx = 7^ is non-zero. 

On the other hand, by the Chern class representa- 
tion we have Q{s) — ci[C)f{s). We further reiterate 
that f{s) is represented by its components fa{s), asso- 
ciated with the g raph l inks a, as described in some de- 
tail in subsection IV A that map a torus = 5^ x 



to M^x ^ a^'^ referred to as 2-dimensional cycles in 
Mx >^ T- The disc A, spanned on s can be viewed as a 
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map from a standard disc = {(a;, y) G 



-y ^ 1} a-nd in the low-temperature limit Eq. (71 1 reproduces 



to M.^x tt^at extends the driving protocol s viewed as a 
map from a circle that represents the disc boundary 
to the space M.'^^ oi robust parameters. That delivers 

X r 



X 



Aix X T that maps 



an extension of the map /q(s) 
to the map faiA) : x - 

X to Aix ^ T- Therefore the cycle /a(s) represents 
the boundary of fa{-A.) [which can be formally stated as 
fa{s) = dfa{A), where d has a meaning of a bound- 
ary operator], i.e., /a(s) can be viewed as the restriction 
of the map fa{A) from its domain x to the do- 
main's boundary x S^. The value of the Chern class 
ci{C) evaluated at a boundary is always zero; this in- 
trinsic property of the Chern classes, which reflects their 
topological (homological) character, is demonstrated for 
our particular case explicitly in appendix [X] This implies 
Q{s) — ci{C)f{s) = ci{£)df{A) = 0, which contradicts 
Q{s) = 7^ 0, obtained earlier. 

D5c c 



We are now in a position to show that Aix 
Aix- Since in subsection 



IV C we have established 



A^x C A^5f' it is sufhcient to demonstrate that x £ D^^ 
with = implies x G -Mx- To achieve that, we 
build explicitly an asymptotically exact, in the /3 — >■ cxi 
limit, approximation for the non-degenerate ground state 
\ux{y- j3)) of the twisted master operator Hxiv, P) for all 
A e T and all y in some small neighborhood of a;. In 
the integer case under consideration our cell is labeled 
uj = (jj'^aa',p) by a pair of distinct nodes, a pair of 
distinct links and an ordering p of the set of links, as 
described in some detail earlier in this subsection. In 
particular in the cell, i.e. for {E,W) E D^, we have 
Ej = Ej, and = Wa' . Let and X'^ be the 
maximal spanning trees that correspond to the degen- 
eracy resolutions Wa < Wa' and Wa' < Wa, respec- 
tively. Then — Qjj'{X^) — Qjj'{X'^), and therefore 
2?t^ = imphes ljj'{X^) — ljj'{X'^), so that both can 
be denoted by ljji{u)). Since the ground state population 
in the neighborhood of x is substantially restricted to 
the the lowest-energy nodes j and /, the existence of a 
unique preferred path Ijji (w) allows the gauge transfor- 
mation trick, developed in subsection |IV C[ to study the 
generating function to be applied. Namely we approxi- 
mate the ground state of the twisted master operator on 
the original graph with its counterpart associated with 
the subgraph represented by the path ljj'(uj). Since the 
aforementioned subgraph is a tree the ground state of the 
associated twisted master operator is represented by a 
gauge transformed Boltzmann distribution, as described 



in subsection IV C which in the low-temperature limit 
is concentrated in the nodes j and /. This results in 
the following explicit expression for the low-temperature 
limit for the (non-degenerate) twisted ground state 

Z^(/?) = e-^^^ +e-^^^', (71) 



the expression for the twisted ground state obtained for 
{E, W) e Mx, i-e., no simultaneous degeneracy, in sub- 
section |IVCJ Finally, we emphasize that we were able to 
obtain the expression of Eq. ( 71 1 due to the existence of 
a unique preferred path Ijji (u) in some neighborhood of 
X, the latter property resulting from 2?^ = 0, and such a 
construction would not work for "D^^ ^ 0. 



Chern class representation of rational quantized 
currents 



In this subsection we present a Chern class representa- 
tion for the pumped current Q(s) for the rational quan- 
tization case. Compared to the integer case, the dis- 
cussed representation for the rational case is much more 
involved. There are two reasons for that. First, Chern 
classes are integer-valued by nature, and to be able to 
reproduce rational numbers some additional framework 
needs to be introduced. Second, in the integer case the 
Chern class representation results from a strong prop- 
erty of the generating function, established in subsec- 
tion |IV C[ the geometric contribution to the generating 
function in the low-temperature adiabatic case becomes 
purely topological, and has a very simple form (as we es- 
tablished explicitly) in terms of the relevant Chern num- 
bers. We analyzed some simple examples with permanent 
degeneracy, and despite the existence of explicit expres- 
sions for average currents, the expressions even for the 
low-temperature adiabatic limit of the generating func- 
tions look complicated. The rational Chern class formula 
for Q(^s), presented in this subsection, should be viewed 
as just a translation of the averaging formula to the lan- 
guage of fiber bundles and their characteristic classes. 
Since in the rational case the discussed representation is 
more involved and less understood, we will restrict our- 
selves to formulating the rational Chern class representa- 
tion and presenting some arguments in its support, post- 
poning the details to some other publication focusing on 
the topological aspects of the stochastic pumping effect. 

An ability to obtain rational numbers is provided by 
bringing in X-theory with rational coefficients [31] . Con- 
ventional topological i^-theory that can be also referred 
to as if-theory with integer coefficients deals with vector 
bundles, i.e., the bundles whose fibers are represented by 
finite-dimensional vector spaces (linear bundles that play 
an important role in our Chern class picture of stochas- 
tic currents represent a particular case of vector bun- 
dles of rank 1). Elements of K{X), associated with a 
topological space Y are vector bundles, or strictly their 
equivalence classes. First of all, isomorphic bundles rep- 
resent the same element of K{Y). There is a coarser 
equivalence relation also imposed, referred to as stable 
equivalence, where we say that two bundles E and E' 
are stable equivalent E ^ E' if there is a bundle F so 
that E (B F = E' (B F , i.e., the bundles become isomor- 
phic after being summed with some bundle F. Note that 
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we can form their direct sum E (B E' and tensor product 
E®E' by applying the operation of a direct sum and ten- 
sor product of vector spaces fiberwise. To convert K{Y) 
to an abelian group, or more precisely a ring we apply 
the Grothendieck group construction by adding to K(Y) 
"negative bundles", so that an element of K{Y) is rep- 
resented by a formal difference E — E' with the stable 
equivalence relation. Equivalence classes [E] — [E'] that 
represent the elements of K{X) are often referred to as 
stable bundles. Note that the construction that involves 
"negative" bundles is analogous to the one used to obtain 
integer numbers Z out of natural ones N. 

Before starting a brief description of the rational K- 
theory we consider a simple example that shades light 
on the importance of the stable equivalence condition. 
Let T(S'^) be the tangent bundle over a 2-dimensional 
sphere, whose fiber TxiS"^) consists of all tangent vectors 
to the sphere that originate at point x £ 5*^, where the 
sphere is assumed to be embedded S"^ C M'^ into the 3- 
dimensional space in a standard way. It is intuitive that 
the tangent bundle Tx{S'^) is non-trivial, i.e., its fibers 
may not be parallelized globally (i.e., for all points of 
the sphere) in a continuous fashion. The normal bundle 
Af{S'^), whose fibers are represented by the vectors or- 
thogonal to the sphere, is trivial though, since it allows 
a global continuous non-zero section, i.e., associating in 
a continuous fashion with any point in the sphere a non- 
zero vector in the fiber, represented, e.g., by a unit nor- 
mal vector directed outside of the sphere. The direct sum 
T{S^)®Af{S^) ^ S^xR^ is, however, trivial, since the to- 
tal space of the bundle is given by a cartesian product of 
the base with the fiber E'^, with the latter represented 
by the 3-dimensional space the sphere is embedded into. 
Therefore, the non-trivial tangent bundle T{S'^) is stable 
equivalent to a trivial bundle. On a more general note in 
the stable equivalence world one can identify the normal 
bundle as Af{Y) = —T{Y) the negative tangent bundle 
for any compact topological space Y, and up to a trivial 
bundle. 

if -theory with rational coefficients associates with any 
space Y an algebra K(Y;Q) over rational numbers, i.e., 
elements of K(Y] Q) can be added, multiplied, and mul- 
tiplied by rational numbers. Elements of K{Y; Q) are 
represented by formal superpositions of complex vector 
bundles over Y with rational coefficients. Similar to the 
integer case, the algebra K{Y;Q) is given by imposing 
the stable equivalence relation on the set of representa- 
tives, which looks absolutely the same as in the integer 
case. Addition and multiplication are naturally defined 
by the operations of direct sum and tensor product of 
vector bundles. Elements of K{Y-, Q) are referred to as 
rational stable bundles. 

Stable equivalence combined with rational coefficients 
allows complete analysis of vector bundles, which means 
rational stable bundles in terms of their Chern classes. 
The Chern classes are powerful topological invariants of 
vector bundles over complex numbers (i.e., the fibers are 
vector spaces over C). However, they generally can- 



not distinguish between any two non-equivalent bun- 
dles. Stable equivalence together with rational coeffi- 
cients turns the situation around, namely a rational sta- 
ble bundle is fully characterized by the set of its Chern 
classes. 

When characterizing the rational stable bundles in 
terms of their Chern classes, it is more natural to 
deal with the Chern characters, rather than Chern 
classes themselves. The Chern character ch{C) of £ e 
K{Y; Q) is a universal combination of Chern classes of 
C, and is represented by a rational superposition of even- 
dimensional cocycles in Y , i.e., linear functionals on cy- 
cles of Y . Stated equivalently a 2fc-cocycle can be eval- 
uated on a 2fc-cycle of Y , resulting in a rational num- 
ber. Note that for example the first Chern class ci is 
a 2-cocycle with integer coefficients. The convenience of 
Chern characters is determined by their algebraic prop- 
erties ch{C ® C) = ch{C) + ch{C') and ch{C ® £') = 
ch{C) ■ ch{C'). For a linear (i.e., 1-dimensional) bundle 
C we have 

ch{L) = e^^(^) = l + c,{C) + ^(ci(£))2 + (72) 

where the the term (ci(£))^ and higher ones should be 
evaluated on the cycles of dimension 4 and higher, re- 
spectively. 

The reason why we brought in if-theory with ratio- 
nal coefficients is that if instead of having one preferred 
line bundle C over A^x x T that describes the pumped 
current Q{s) via its first Chern class ci(£), as is in the 
integer case, we have a set of line bundles, and we can 
average them with rational weights, resulting in an "av- 
eraged bundle" C G K{M.x x 7";Q) and further express 
the pumped current via its Chern character ch{C). 

Recalling the definition of 2-cycles fa{s) in the space 
A4x xT, labeled by the links a G Xi of our network X, 
the Chern class representation of Q{s) has the form 



Qa{s) ^ ch{C)fa{s), C= y2 ^kX' 

nx 

{k,X) 



(73) 



where C E K{Mx x T; Q), and a construction of the line 
bundle £^ , associated with a degeneracy resolution, 

labeled by {k,X) is presented in appendix [b] 

The Chern class representation [Eq. ( 73 )| can be jus- 



tified by recasting it in an equivalent form 
1 

nx 



E 



Cl 



C 



k,X 



fa{s), 



(74) 



followed by applying the definitions and the averaging 
formula. 

The extension of the Chern class representation 



Eq. (73)] to the extended subspace Mx 3 Mx is as 



follows. The above inclusion generates an algebra homo- 
morphism K{Mx x T;Q) ^ K{Mx x T;Q) that boils 
down to restriction of a bundle to a subspace. The space 
satisfies the following properties: (i) the rational 
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stable bundle £ G K{Mx x T;Q) allows a unique ex- 
tension to a stable bundle € iir(A^5f x T; Q), and (ii) 
cannot be extended to a rational stable bundle over 
U xT for any open U D Aix and U ^ A^x- The space 
Mx is fully characterized by above properties, and so is 
C' . For a driving protocol s that stays entirely in 
the Chern class representation reads 

Oa(s) = c/i(£'^)/„(s). (75) 

It is also instructive to note that since the line bun- 
dles £^ are trivial over Ue and Uw , so is there direct 
sum nxC = ®(fcx)''^fcx- This means that nxC can 
be viewed as a complex vector bundle of dimension nx 
obtained by gluing together the trivial bundles defined 
over Ue and Uw, so that the Chern class representation 
for the current can be recast in a form 

Qc.{s) = —ciinxC)Us). (76) 
nx 

VI. DISCUSSION 

We have demonstrated that topological invariants, 
such as Chern classes/numbers, naturally appear in clas- 
sical stochastic processes when we consider counting 
statistics of generated currents under the condition of pe- 
riodic time-dependence of kinetic rates. We determined 
that the integer quantization of currents, observed in the 
adiabatic low temperature limit, is a special manifesta- 
tion of such Chern numbers. The regime of fractional 
quantization remains much less understood. We calcu- 
lated the counting statistics in this regime for a simple 
model with 2 nodes and N links. The result appears rem- 
iniscent the averaging formula, in which a weighted sum 
of the reference integer circulating (conserved) currents 
Q, associated with legitimate degeneracy resolutions, as 
described in detail in the first manuscript of the series, 
is replaced by the weighted sum of the e*-^*^ terms with 
the same weights. We leave the question of possible gen- 
eralization to the case of arbitrary graphs as an open 
problem. 

We also demonstrated that the space of robust pa- 
rameters for cyclic protocols that generate the integer- 
quantized currents, which are robust with respect to 
small perturbations of the driving protocol, is the same 
as the space of parameters that allow the standard adia- 
batic approach to be applied to the generating function 
calculation for all values of its argument. This results 
in the identification of a topological invariant (based on 
first Chern class) of the counting statistics that is directly 
related to the quantized average current. The described 
relation provides a simple approach to classifying topo- 
logical properties of the counting statistics. 

Our results are only a first step towards understand- 
ing the topologically protected behavior of the counting 
statistics. We focused on the low-temperature adiabatic 
Q{s) = lim^_>oo limT-^_j.oo 0('S; 'T'd; /?) limit. Neverthe- 
less, we showed that topological invariants in counting 



statistics can be defined for arbitrary temperatures and 
frequencies. In this more general situation they no longer 
directly lead to the quantization of average currents, but 
still may correspond to some other effects. 

Other types of topological invariants in counting statis- 
tics, which have not been discussed in the present article, 
are possible. For example, in [21 fractional quantiza- 
tion of currents in a model without detailed balance was 
found. When detailed balance is not in place the twisted 
master operator is no longer Hermitian. The correspond- 
ing geometric phases may have different topological prop- 
erties due to existence of the so-called exceptional points 
|13j . Other topological invariants in counting statistics, 
called Z2-invariants, have been recently found and stud- 
ied [ig. 

We would also like to note that a direct relation 
between the twisted master operators and Hermitian 
Hamiltonians suggests that all of the phenomena we have 
discussed in this series in terms of applications to stochas- 
tic processes can be reinterpreted in terms of the cor- 
responding properties of quantum mechanical systems. 
For example, in the present manuscript we have provided 
a simple and intuitive procedure for calculating integer 
topological invariants of the operator eigenstates phase 
space, parametrized by the periodic variables A and s. 
In quantum mechanical applications, such topological in- 
variants may describe Chern classes of vector bundles, 
associated with Bloch bands. Hence, the simple classi- 
fication approach, based on the averaging formula, may 
appear useful for studying topological properties of some 
Bloch bands. 

In both manuscripts of the series we have been study- 
ing generation of currents for stochastic motion in net- 
works, or in other words on graphs, i.e., in the discrete 
setting. Similar questions can be asked for continuous 
setting, i.e, in the case of Langevin stochastic processes 
in continuous spaces, e.g., manifolds, where the role of 
the master equation/operator is played by its Fokker- 
Planck counterpart. When the configuration space X 
where stochastic dynamics occurs has non-contractible 
1-dimensional cycles (loops) global stochastic currents 
(fluxes) can be introduced (see, e.g., [35]) and observed in 
experiment. The detailed balance condition corresponds 
to potential forces that describe effects of advection, and 
driving can be described by time dependence of a po- 
tential function V{x). A generic potential function is 
of Morse type, i.e., its critical points are isolated and 
have non-degenerate matrices of second derivatives. In 
the low-temperature /3 — > c» limit, the Boltzmann distri- 
bution p^{x) = Z~^e~^^^^' is concentrated in a small 
neighborhood (the lower the temperature, the smaller the 
neighborhood) of the global minimum of the potential 
V{x). There are two mechanisms that generate current 
in the continuous case. In the first case, consider a situa- 
tion when the potential evolves in time in a way that the 
positions Xj of the local minima stay the same, while the 
local minimal values V{xj) of the potential function are 
subject to changes in time. When the global minimum 
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FIG. 5. Contour of integration 



at, say Xi, loses its global character, delegating it, say 
to ajfc, a current will be produced due to the probability 
flux from a neighborhood of Xi to the neighborhood of 
Xk, with the flux concentrated around the minimal en- 
ergy path, the latter generally passing between a set of 
some other local minima via elementary paths that go 
through the transition states, i.e., critical points ya with 
a single unstable mode. The described process can be 
readily analyzed by applying the approaches, developed 
in this series, to a network, whose nodes and links are 
represented by Xj and ya, respectively, while the rates 
are computed using standard transition rate theory for 
overdamped Langevin dynamics (see, e.g., [40] ) . The sec- 
ond mechanism for generating current is just simple pe- 
riodic motion of a global minimum. Strictly speaking, 
current generation is a combination of the above two 
mechanisms. 

Most of the structures of our present study can be 
extended to the continuous case using more technical, 
rather than conceptual effort. In particular, a space of 
good parameters that has a property of describing generic 
potential functions can be relatively easily introduced. 
The question of identifying the robust parameter space, 
in the sense of our present consideration, is more com- 
plex. Current quantization for Langevin processes in con- 
tinuous spaces, which is currently under study, will be 
presented in our further publications. 

Finally, the fact that certain properties of fluctuations 
in mesoscopic systems can be robustly controlled by pe- 
riodic driving protocols should be of interest for charac- 
terization and control of mesoscopic structures, such as 
electronic nanoscale conductors 136 and molecular mo- 
tors EH El] . 



When computing a generating function, we deal with a 
family of twisted master operators, Hx{x), where A G T 
is the argument of the generating function, whereas 
X represents a set of parameters, e.g., the node ener- 
gies Ej and the barriers Wa, i.e., x = {E,W) with 
E = {Ej\j £ Xq) and W = {Wa\a e Xi). Note that the 
operator Hj^{x) becomes Hermitian upon the transfor- 
mation, Hx{x) — >■ e^^/'^Hxlx)e'^^/'^ , however, its com- 
ponents are generally not real. This means that eigen- 
values of Hx{x) are real, whereas the eigenvectors have 
generally complex components. An adiabatic calculation 
of the generating function for a given driving protocol 
is straightforward, as long as the ground state of Hx{x) 
is non-degenerate for all A g T and all values x of pa- 
rameters involved in the driving protocol [lOj [13] . We 
reiterate that we refer to such parameters sets as non- 
degenerate ground-state parameters. The (temperature- 
dependent) space of such parameters x = [E, W) has 
been denoted by M'^xiP)- Actually, one of the main re- 
sults of this manuscript is the demonstration that in the 
low-temperature adiabatic limit the average current is 
quantized and is robust with respect to perturbations of 
the driving protocol as long as the latter is restricted to 
the subspace, A^5f ^ M.x, of robust parameters, which 
is a kind of a low-temperature limit space for Ai'^xil^)'^ 
the precise construction is presented in subsection |IV B| 

Having restricted ourselves to the subspace of non- 
degenerate ground-state parameters, we note that with 
each point y ~ {x, A) = [E, W; A) we can associate 
a complex line (complex 1-dimensional vector space) of 
eigenstates such that Hx{x)\u) = Wyj?/), where ujy 
is the largest eigenstate of Hx{x). Considering the de- 
scribed above eigenstates for all Hx{x) we obtain a space 
£, referred to as the total space of a vector bundle, that 
consists of points (a;, A, \u)). The space of eigenstates, 
\u), can be viewed as fibers of the vector bundle, and are 
parametrized by {x,x)- 

We are now in a position to introduce the flrst Chern 
class ci(£) of £. We start with noting that the fiber Cy 
over y = {x, A) consists of the states \u) that differ by 
just normalization factors, and for any y we can locally, 
i.e., in some neighborhood of y, fix the normalization 
factors in a smooth way, i.e., build a smooth function 
\ux{x)) of {x,\) with Hx{x)\ux{x)) = ujy\ux{x)). Such 
a function is referred to as a local section of C. We further 
introduce a vector-potential 



Appendix A: First Chern class and relevant Chern 
numbers 



My) = 



{u{y)\dju{y)) 
{u{y)\u{y)) 



dj = d/dy 



'3' 



(Al) 



In this appendix we present a simple and self-contained 
description of the first Chern class and provide all nec- 
essary information for our applications. Chern classes Cj 
with j = 1,2, .. . are useful invariants of complex vec- 
tor bundles. In our applications we will be dealing with 
linear, i.e., dimension one bundles, which are fully char- 
acterized by their first Chern class ci. 



where we introduce the left eigenstates, defined by 
{u\Hx{x) = {u\ujy, and the scalar product is defined as 
the sum of the product of corresponding components of 
the left and right eigenstates. We define the curvature of 
the vector-potential by 



F,,iy)^d,A,{y)-d,A,{y), 



(A2) 
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the latter satisfying a property 

^^FJk{y) + dkF„{y) + djFk,{y) = 0. 



(A3) 



Eq. (A3 1 is usually referred to as the property of Fij{y) 
to be closed. It is very important to note that the vector- 
potential depends on the particular choice of a local sec- 
tion, i.e. for = e"^^^^ we have 



Ar{y)^A,iy)+dMy), 



(A4) 



however, Eq. (A4) represents a gauge transformation, so 
that the curvature Fij{y) is not sensitive to the local 
section choice, and there fore Fij{y) is defined globally 
and depends only on £ and the scalar product. 

Being interested in applications that involve integra- 
tion of the curvatures over 2-dimensional cycles we define 
that two closed forms, Fij{y) and F-j{y), are equivalent 
if Pljiy) = Ftjiy) + diajiy) - dja^{y) for some glob- 
ally defined smooth ai(y), and hereafter use the notation 
[F] for the equivalence class of Fij{y). The equivalence 
class [F] does not depend on a particular choice of the 
scalar product, so that the first Chern class, defined by 
ci(£) = [F], depends on C only. We will also demon- 
strate that the first Chern class is the same for equiva- 
lent bundles £' ~ £, e.g., when the fibers of £' can be 
continuously deformed to the fibers of £, so that the first 
Chern class Ci{C) is a topological invariant of C. 

We conclude this brief description of the first Chern 
class of line bundles with showing how to evaluate the 
Chern class on 2-dimensional cycles in the parame- 
ter space represented by points y — (a;. A), resulting 
in Chern numbers, followed by demonstrating that all 
Chern numbers are represented by integers. By anal- 
ogy with "standard" 1-dimensional cycles, we define a 
2-dimensional cycle as a compact oriented 2-dimensional 
manifold mapped in a smooth way to the parameter space 
of y. We will restrict ourselves to the case when the do- 
main manifold is a 2-dimensional torus with coordi- 
nates < T < 1 and < X < 27r, also denoted by (C^, C^), 
so that a cycle / is represented by a smooth doubly pe- 
riodic map y(^). The Chern number ci{C)f, which can 
be viewed as a result of evaluation of the Chern class on 
the cycle /, is given by 



ci(/:)/= / dT 



2-iT 



27ri 



where e^^ 



^21 



1 is a 2-dimensional antisymmetric 



object, and hereafter we assume summation over repeat- 
ing indices. For our applications, r will represent the 
dimensionless time on the driving protocol and x is a 
phase of some component of a multi- variant argument of 
the generating function. 

To establish the integer nature of Chern numbers 
we pull-back the bundle C to along / that re- 
sults in a bundle over T^, denoted /*£, whose fiber 



Following the general procedure de- 



scribed earlier [Eqs. (All and (A2|] we construct a 



gauge field Aa{^) and curvature Fab{^). Due to covari- 
ance we have Aai^) = {dy^ /d(a)A/{y{^)) and Fabi^) = 
{df/d^a){dv'/d£,b)F,,{ym^ so that 



ci(/:)/= / dT 



2m 



Fabii)- 



(A6) 



We further view a torus as a rectangle [0, 1] x [0, 2tt] 
with the opposite sides being identified, and build a sec- 
tion of f*C over the rectangle by fixing the normalization 
of \u{0,x)) in some arbitrary way (which is obviously al- 
ways possible to do), followed by extending the normal- 
ization to the whole rectangle by imposing the condition 
(m(t, x)|9t-u(t, ^)) = 0. This results in a section |'u(^)), 
which is quasi-periodic, i.e., |u(r, x-|-27r)) = |u(t, x)) and 

|u(t + 1,x)) =C(x)|u(t,x)), withC(x + 27r) = C(x)- Due 
to the quasi-periodic nature of |u(^)) the vector poten- 
tials on the opposite sides of the rectangle differ by just 
gauge transformations [Eq. (A4)]: 



Ar{T,2n) = Ar{T,0), 

4(l,x)-4(0,x) + a;,C(x)/C(x) (A7) 

Integrating the vector potential over the border of the 
rectangle in the counterclockwise direction we arrive at 



1 

27h 



d^AM) 



dx dCix) 
27r<(x) dx 



n[ax)m 



where ?i[C(x)] is the number of times zrj C winds over 
the zero point in the complex plane C while x varies from 
to 27r. Since Stokes' theorem (see e.g., !42 ) identifies 
the l.h.s. of Eq. (|A8|) with the r.h.s. of Eq. (|A6|), we can 



express the Chern number as 



ci(/:)/ = n[c(x)]. 



(A9) 



Where the r.h.s. of the later is the winding index. We 
wish to emphasize that, from its construction, C(x) raay 
be interpreted as the adiabatic geometric factor associ- 
ated with a cycle of adiabatic dynamics with a twisted 
master operator H\, where A = e*-^. Therefore, C(x) 
describes the dependence of the geometric factor on the 
argument of the generating function, whereas the Chern 
number ci{C)f represents an integer topological invari- 
ant of C(x)j and Eq. (A9| will play a key role in identify- 
ing the Chern numbers as the components of the quan- 
tized currents in the low-temperature adiabatic regime of 
stochastic pumping. 

In subsection |V A| we have used the intrinsic property 
of Chern classes to vanish being evaluated at boundaries. 
Here we demonstrate this property explicitly for our par- 
ticular application. For a map A : — > where 
i s a s tandard 2-dimensional disc defined in subsec- 



tion 



VA 



we have a map fa 



X r. 



defined by fa{A{r],x)) = (-^(77), A" (x)), where the com- 



ponents of A"(x) are: A2(x) 



and A"(x)= 1 



for 7 ^ a, as described at the end subsection |IV A| in 
the contest of defining the set / of 2-dimensional cycles. 
The map faiA) is defined in a way so that the cycle 
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fa{s) = dfa{A) represents its boundary, i.e., fa is ob- 
tained from fa{A) by restricting the domain x to 
its boundary d{D^ x S^) — x S^. We further apply 
the expression of Eq. (|A6|) 



Ci{C)fa{s) 



d^^e-'Fabi^), (AlO) 



written in terms of some coordinates of the torus. We 
further note that the curvature is actually defined over 
X ^1 by FabiC) = {drJdCa){dR/dCb)FMi)), and 
being restricted t o its b oundary x reproduces the 
curvature in Eq. (AlO). This allows the general (multi- 
dimensional) Stokes formula to be applied 



F= dF, 



(All) 



since x = d{D'^ x S^). The Stokes formula, pre- 
sented in Eq. (All) using the language of differential 
forms ([12] [15]) has the following coordinate represen- 
tation 



Six Si 

3 /• ^abc 

u 

D^xSi 



dXe'^'^dahciO. 



where e is the 3-dimensional antisymmetric object. 



Due to the definition of the curvature [Eq. (A2)] the cur- 
vature is closed, i.e., dF — 0, or in coordinates is given 
by Eq. ([A SJ). Therefore the integrand in the r.h.s. of 



Eq. (A12) turns to zero, which combined with Eq. (AlO) 



yields ci{C)dfa{A) — 0, i.e., vanishing of the Chern class 
evaluated at a boundary cycle. 



Appendix B: Construction of the current generating 
stable bundle 



According to Eq. ([73]) the current generating stable 



allowed degeneracy resolutions. In this appendix we 
present a construction of the line bundle Cj^ j associ- 
ated with a degeneracy resolution, labeled by (fc,X). 
The construction is based on gluing together a set of 
bundles defined over the open sets that cover the bundle 
base Mx x T . 

The subspace A^x is obviously covered by a family of 
open sets 



U y C/[vv], (Bl) 



where U\e\ and consist of all {E\ W') so that 

\E\ = \E\ and [VK]' = \W\ respectively. Note that 

v/n = for \E\ ^ [£;]' 



\E'\ 



= and Uyw\ n U\ 



and \W\ 7^ [W']', respectively. In particular there are 
no triple or higher-order intersections among the sets 
of the family. Given a degeneracy resolution ^fc, X), 
(A12) '^6 can build complex line bundles c}-^^ and ' over 

\E\ 



bundle £ is obtained via averaging of £. js- over all 



V{E\ X T and C/[vi/] x T, respectively. The fiber of 
over [E, W, A) G U[e] x T is spanned on the state Ifc^;), 

whereas the fiber of -C^^' over {x, A) = (E, W, A) S 
U[w] X T is represented by the eigenstates \u) of the op- 
erator H\{x; P)\x^ with the lowest eigenvalue. The op- 
erator Hx{x; ji) I is obtained by restricting the twisted 

master operator H\{x\ /3) to the spanning tree Xw^ i.e., 
the rates associated with the withdrawn links are set 
to zero. The fibers of C^^^ and ^C^^' over all points 
(a;. A) e {U[E] n U[w]) x T should be glued together. 
This is done by just projecting an eigenstate \u) onto 
jfe^;). Since there are no triple intersections of the open 
subsets, no consistency conditions have to be taken care 
of. 
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